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ABSTRACT 


In  the  absence  of  wave  damping  the  momentum  of  a weak  cold  beam  can 
be  clamped  while  the  wave  amplitude  increases  secularly  in  time.  When 


wave  damping  is  present  a nonlinear  equilibrium  is  attained  in  which  the 


beam  evolves  into  a singular  charge  clump  which  drifts  at  constant 


velocity  through  the  plasma 


k 


It  is  well  known  that  when  a low  density  fast  electron  beam  is  injected 
into  a background  plasma  it  triggers  an  electrostatic  two-stream  instability. 
It  has  been  shown  both  theoretically1’2  and  experimentally3’4  that  this 
instability  leads  to  the  growth  of  a nearly  monochromatic  electron  plasma 
wave  (single  mode)  whose  amplitude  saturates  by  trapping  the  beam  electrons 
within  its  potential  troughs.  The  process  of  beam  trapping  is  of  interest  to 
basic  plasma  physics  studies  since  it  provides  an  example  where  coherent 
nonlinear  processes  play  a dominant  role,  thus  entering  the  domain  where 
the  familiar  quasilinear  theory  based  on  the  random  phase  approximation  is 
not  applicable.  On  the  practical  side,  saturation  by  beam  trapping  limits 
the  maximum  energy  that  can  be  extracteu  in  travelling  wave  tube  (TWT) 
devices  used  extensively  to  generate  microwave  radiation. 

This  letter  illustrates  the  changes  produced  on  the  dynamics  of  beam 
trapping  when  an  external  DC  electric  field  Eq  is  applied  to  the  beam 
particles.  This  problem  is  of  interest  for  various  reasons:  (1)  it  provides 
an  idealization  of  the  coherent  interaction  of  a population  of  runaway  elec- 
trons with  collective  modes,  (2)  it  explores  the  possibility  of  enhancing 
the  efficiency  of  energy  extraction  out  of  beam  devices,  (3)  it  simulates 
the  effect  of  low  frequency  turbulence  on  beam  nonlinearities.  The  major 
effects  found  in  this  study  are:  (1)  the  beam  can  be  clamped  in  energy 
with  the  momentum  push  transferred  to  the  wave,  (2)  the  normal  beam  satur- 
ation level  can  be  significantly  enhanced  by  Eq  ^ 0 (DC  to  AC  conversion), 

(3)  a dynamic  Bemstein-Greene-Kruskal  (BGK)  equilibrium5  is  attained  when 
wave  damping  is  included. 

The  present  study  is  based  on  the  formulation  of  the  small  cold  beam 
problem  given  by  O'Neil,  et  al.1  The  essence  of  the  procedure  consists  of 
describing  the  behavior  of  the  background  plasma  analytically  through  a 


linear  dielectric,  and  following  the  exact  orbits  of  the  beam  electrons 
in  a computer,  i.e.,  the  beam  is  discretized  into  N particles.  The  velocity 
Vj  and  position  of  the  j particle  is  obtained  from 

” 5t  vj  3 eEo  " eE(t)exp{i[kXj(t)  - a)Qt]  > + c.c.  (1) 

where  e,  m refer  to  the  electron  charge  and  mass,  k is  the  wavenumber  of 
the  fastest  growing  wave,  tuQ  is  the  frequency  of  the  wave  in  the  absence  of 
the  beam,  and  E(t)  = E(0)exp[-  i /^dt * 6o>Ct * ) ] refers  to  the  slowly  changing 
wave  amplitude  with  being  the  complex  frequency  shift.  The  wave  amplitude 
is  determined  self-consistently  from  Poisson's  Eq. 


i)  Ife  ■ -mr  ‘fyt)  - v»  (2> 

' o J 

where  n^  is  the  beam  density  and  e the  linear  dielectric  of  the  background 
particles,  i.e.,  e(k,  o)q)  = 0. 

The  physics  associated  with  Eq  t 0 is  made  clear  by  the  two  conservation 
laws  associated  with  Eqs.(l)  and  (2) 
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with  Uj  = Vj  - wQ/k.  For  Eq  * 0,  Eqs.(3)  and  (4)  are  satisfied  by  balancing 
the  increase  in  wave  amplitude  with  the  simultaneous  recoil  and  trapping  of 
the  beam.  However,  for  Eq  ^ 0 Eq.  (3)  indicates  that  the  total  momentum  of 
the  beam-wave  system  must  increase  secularly  in  time.  There  are  two  ways 
to  satisfy  this  constraint:  (1)  the  beam  momentum  increases  secularly  in 
time  with  the  wave  amplitude  essentially  constant,  i.e.,  the  beam  runs  away. 
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and  (2)  the  bean  momentum  remains  clamped  while  the  wave  amplitude  increases 
secularly.  A given  state,  beam  or  wave  runaway,  is  determined  by  the  value 
of  Eq.  The  threshold  value  separating  the  two  states  is  obtained  naturally 
from  this  formalism  because  it  can  be  fully  scaled.  After  going  through  the 
scaling  the  equation  of  motion  shows  that  beam  clamping  is  possible  if 
Eo  < ET*  n 3 [(nb/n0)  “J/C^/3u,)o]l/3,  and  is  the  electron 

plasma  frequency  of  the  background  plasma  of  density  nQ.  Ej,  is  the  saturation 
amplitude  for  Eq  = 0.  Physically,  if  Eq  > Ej.  the  wave  can  not  cancel  the 
external  force,  thus  leading  to  a runaway  beam.  For  Eq  < Ej,  the  external 
push  prevents  the  beam  from  recoiling,  thus  the  wave  continues  to  be  driven 
by  the  beam. 

To  verify  the  existence  of  the  two  runaway  states  and  their  dependence 
on  we  have  solved  numerically  for  the  time  evolution  of  the  system 

using  the  scaled  versions  of  Eqs.(l),  (3),  (4).  Figure  1 exhibits  the  time 
evolution  of  the  scaled  amplitude  (A  = E/ET,  t = ftt)  and  the  momentum  of  the 
beam  P = £ ku^/Ntt.  The  DC  field  is  turned  on  at  t = 10.0  at  a value 
E /Ej.  = 2.0.  It  is  seen  that  the  amplitude  saturates  at  a level  slightly 
larger  than  that  for  Eq  * 0,  but  it  does  not  execute  the  typical  trapped 
particle  oscillations  because  all  the  beam  electrons  are  untrapped,  as 
indicated  by  the  secular  increase  in  P.  In  this  case  the  beam  runs  away, 
as  expected,  and  the  wave-particle  resonance  is  detuned  for  this  particular 
mode.  Of  course,  other  modes  with  smaller  k are  expected  to  be  excited  as 
the  beam  accelerates. 

Figure  2(a)  shows  the  time  evolution  of  the  wave  amplitude  for  the 
case  Eo/Ej.  = 0.5.  As  predicted,  the  wave  is  observed  to  grow  secularly  to 
large  amplitude.  The  Eq  = 0 result  is  included  for  comparison.  Figure  2(b) 
shows  the  evolution  of  the  frequency  shift  W.  An 

■ 
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examination  of  phase-space  (not  shown)  indicates  that  the  bulk  of  the  beam 
electrons  (>  90%)  is  trapped  within  the  troughs  of  the  growing  wave,  as  is 
evident  in  the  amplitude  oscillations  of  Fig.  2(a).  In  addition,  a small 
secondary  runaway  beam  is  created;  it  arises  from  the  few  particles  (<  10%) 
which  spill  over  the  potential  well  as  the  beam  rotates  in  phase-space. 

As  indicated  by  the  conservation  laws,  the  beam-wave  system  can  not 
attain  a steady  state  for  Eq  f 0.  To  reach  a steady  state  an  external 
dissipation  mechanism  must  be  present,  e.g.,  plasma  heating,  wall  heating, 
radiation  leakage.  To  investigate  this  effect  we  introduce  a scaled 
damping  factor  T/2  whose  origin  need  not  be  specified.  The  scaled  conser- 
vation laws  become 

ft  [P  + A2]  = (Eo/Et)  - TA2  (5) 

^ [K  ♦ 2WA2]  - (Eq/Et)P  - 2WTA2  (6) 

where  W is  the  real  part  of  the  scaled  6u  and  K is  the  scaled  kinetic  energy 
of  the  beam.  A steady  state  solution  of  Eqs.(5)  and  (6)  is  A = (E^/E^r)55, 

P = W = 0.  However,  to  satisfy  P ■ 0 for  all  t,  the  equation  of  motion 
demands  that  EQ/ET  = A sin  (kx^-Q)  for  some  fixed  phase  factor  6.  Since 
this  constraint  must  hold  for  all  j it  follows  that  all  beam  particles  must 
occupy  the  same  position,  i.e.,  the  beam  must  evolve  into  a singular  charge 
clump  whose  momentum  remains  fixed  although  an  external  electric  field  is 
being  applied.  This  unusual  nonlinear  state  is  an  example  of  a dynamic 

BGK  equilibrium  in  which  the  beam-clump  acts  as  an  intermediate  agent  that 

transfers  the  momentum  push  of  Eq,  through  the  wave,  to  those  mechanisms 
that  give  rise  to  r.  The  self-consistency  of  the  clump  position  and  the 
amplitude  demand  that  (EQr/E^)  < 1 for  this  steady-state  to  exist. 


1 
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Figure  3 illustrates  the  evolution  of  A and  W for  a case  Eo/ET  * 0.5, 
r = 0.05,  satisfying  the  existence  criterion.  It  is  seen  that  the  runaway 
wave  saturates  and  approaches  asymptotically  the  predicted  steady  state 
value  of  •'To.  The  frequency  shift  approaches  W ■ 0,  as  is  required  for  the 
steady  state.  The  evolution  of  the  clamped  beam  toward  a singular  charge 
clump  is  demonstrated  in  the  phase  space  plots  shown  in  Fig.  4.  For  t = 25.0 
we  see  the  pattern  associated  with  the  early  wave  runaway  process  typical 
of  the  r = 0 behavior  shown  in  Fig.  2.  However,  for  late  times  (x  = 85.0) 
the  bulk  of  the  beam  particles  form  a singular  clump,  as  predicted.  It  is 
also  observed  that  a small  secondary  runaway  beam  is  generated,  but  its 
presence  does  not  affect  the  attainment  of  the  singular  equilibrium.  Of 
course,  physically  one  expects  that  the  singular  charge  clump  is  limited 
in  extent  in  a real  system  due  to  the  strong  repulsive  forces  that  arise, 
such  an  effect  is  not  included  in  the  present  formulation. 

In  summary,  the  present  study  shows  that  when  the  exact  beam  dynamics 
are  followed  it  is  possible  for  a runaway  beam  to  become  clamped.  This 
result  confirms  an  earlier  prediction6  of  this  effect  based  on  a spatially 
averaged  formalism.  The  clamping  process  generates  a runaway  wave,  thus 
suggesting  an  efficient  method  for  converting  DC  to  AC  energy.  When  wave 
damping  is  present,  a dynamic  BGK  equilibrium  is  approached  in  which  the 
beam  evolves  into  a singular  charge  clump  that  drifts  through  the  plasma 
at  constant  velocity.  The  stability  and  limiting  factors  of  this  state 
require  further  study. 

The  author  thanks  Dr.  G.  Dimonte  for  useful  conversations  on  this  topic. 

Work  supported  by  the  Office  of  Naval  Research  and  USDOE  Contract 
No.  EY-76-C-03-0010  PA  26. 
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FIGURE  CAPTIONS 


Time  evolution  of  the  scaled  wave  amplitude  A and  average 
beam  momentum  P in  the  beam  runaway  regime  (Eq  > E^,) , for 
Eo  = 2.0  EJ  , r = 0. 

(a)  Time  evolution  of  the  scaled  wave  amplitude  A and 

(b)  frequency  shift  W in  the  clamped  beam  regime  (Eq  < ET) , 

E = 0.5  E.  , T = 0.  Result  for  E = 0 is  included  for 
o T o 

comparison. 

Evolution  toward  the  dynamic  BGK  equilibrium  in  the  presence 
of  wave  damping,  EQ  = 0.5  , r = 0.05. 

Phase-space  plot  showing  the  formation  of  a singular  charge 
clump,  Eq  = 0.5  , T = 0.05  , (a)  t = 25.0  , (b)  t = 85.0. 
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